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Abstract 

The bipolar nebula, IRAS 19312+1950, is a unique SiO maser source exhibiting both properties as 
young and evolved objects. To clarify the nature of this object, we made molecular line observations with 
the Nobeyama 45-m radio telescope. We detected emission from O-bearing (HCO + , SiO, SO, and SO2); 
C- and N-bearing molecules (CN, CS, HCN, HNC, NH 3 , N 2 H+, HC 3 N, H 2 CS, and CH 3 OH), and their 
isotopic species (C 17 0, 13 C 18 0, and C 34 S). Line profiles consist of a weak broad (Av ~ 30 km s _1 ) and/or 
a strong narrow (At; < 5 km s _1 ) components, depending on molecular species. Strong time variations of 
H2O and SiO masers were also observed. Numerical modeling of the envelope with the LVG-code results 
in a good fit of the model with a mass loss rate of 2.6 x 10 _4 M Q yr _1 to the observed intensities for the 
broad-component lines. On the other hand, non-O-bearing molecules, which have only the narrow profiles, 
were found to have abundances typical to those in cool dust clouds. No isotopic enrichment was found, 
indicating little evidence of the narrow cool component being an ejecta of the central AGB star or a possible 
companion. These facts compelled us to conclude that IRAS 19312+1950 is an exotic mass-losing evolved 
star embedded in a low-mass (~ 20 Mq) dark cloud. 
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1. Introduction 

IRAS 19312+1950 exhibits a prominent bipolar neb- 
ulosity of the size of about 30" in near-infrared im- 
ages (Nakashima, Deguchi 2000). The IRAS color of 
this source, log(F 2 5/ ' F 12 ) = 0.47, indicates a presence 
of cool dust of a temperature of about 200 K toward 
this source. Detecting SiO and H2O maser lines in 
this object, Nakashima, Deguchi (2000) concluded that 
this is an oxygen-rich protoplanetary nebula similar to 
OH231.8+4.2. The distance to this object was estimated 
to be about 2.5 kpc (Nakashima et al. 2004; hereafter 
Paper I), based on the near- and mid-IR flux densities. It 
is slightly unusual that CO and the other thermal molec- 
ular lines of less abundant species are detectable in the 
envelope of an Asymptotic Giant Branch (AGB) or post- 
AGB stars at 2.5 kpc away from the sun. For example, the 
distance to the similar object, OH231.8+4.2, is about 1.3 
kpc (Morris et al. 1987), and those of other commonly ob- 
served AGB stars in radio are less than 1 kpc (Olofsson et 
al. 1998). Near-IR imaging of this object with the UH 88- 
inch telescope shows a bright central star and surrounding 
nebula, extending from NE to SW by about 30" (figure 
1) . It also shows a small (~ 5") ring- like structure near 
the central star, which is elongated perpendicularly to the 



bipolar axis. 

In paper I, we reported the results of mapping obser- 
vations by several molecular lines (CO, HCN, and CS) 
using the Nobeyama 45-m radio telescope; the line pro- 
files exhibit broad and/or narrow components depending 
on molecular species. We found that both broad and nar- 
row components centered at Vi sr ~ 35 km s^ 1 are strongly 
peaked in intensity at the IRAS position of this source 
with a 16" telescope beam. Without doubt, the broad 
component with a full width of more than 20 km s _1 must 
originate from the outflowing envelope of the central star. 
However, the spatial association of the narrow component 
to this object is a puzzle. The mass of the narrow compo- 
nent was estimated to be about 10 M@ (within 16" beam). 
The mass obtained is compatible with the core masses of 
dark clouds (for example, Fuller et al. 1991). 

In this paper, we present a part of molecular-line obser- 
vations toward this object, and discuss on it from a view 
point of circumstellar chemistry. In usual cases, molecu- 
lar species found in dark clouds are considerably different 
from those found in the evolved objects [see, Glassgold 
(1996) and Ehrenfreund, Charnlcy (2000)]. Therefore, the 
origin of the cloud surrounding the central star of IRAS 
19312+1950 can be resolved from this study in some de- 
gree. Although evidence of an evolved star is still not 
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Fig. 1. JH K false-color composite image of the field of IRAS 
19312+1950, using the UH 88-inch telescope with the SIRIUS 
near-IR array camera (by courtesy of M. Tamura and K. 
Murakawa, National Astronomical Observatory). The image 
size is 128" X 128". The north is up, and the east is left. 

completely crucial for this object or rather conflicting, we 
believe that this is a unique evolved object surrounded by 
a remnant of raw material. We present the observational 
results and discussion on the nature of this object in this 
paper. 

2. Observations and results 

The observations were made with the 45-m telescope at 
Nobcyama during 2001 February 4-9, and during 2002 
December 11-21. A few additional observations were 
made on 2001 March 19 and 2002 May 19 for checking time 
variation of SiO maser line intensities. A cooled SIS-mixer 
receiver (SlOO) with a bandwidth of about 0.5 GHz were 
used, and the system temperature (including atmospheric 
noise) was 300-500 K (SSB) depending on the frequency 
and weather. In addition, HEMT and SIS-mixer receivers 
(H22 and S40) were used for the 22 GHz H 2 and the 
43/49 GHz SiO/CS line observations, respectively. The 
half-power beam widths, and the aperture and beam effi- 
ciencies of the telescope at the observed frequencies were 
summarized in table 1. The antenna temperature (T*) 
given in the present paper is that corrected for the atmo- 
spheric and telescope ohmic loss but not for the beam nor 
aperture efficiency. Observations were made in a position 
switching mode, and the off-position was chosen 12' away 
from the object in right ascension, where contamination 
by background CO emission was minimized. Further de- 
tails of the observations were described in Paper I. 

As described in Paper I, strong CO emission was seen 
toward this object. The 15"-, 30"-, and 60"-grid map- 
pings revealed that emission spikes at Vi sr = 27-33, and 
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H22 22 


72 


0.81 


0.62 


S40 43 


40 


0.77 


0.57 


SlOO 86 


18.6 


0.49 


0.42 


SlOO 100 


16.3 


0.51 


0.41 


SlOO 110 


15.4 


0.48 


0.35 



43 km s _1 , are extended more than a few arcmin, in- 
dicating that they are from fore/background clouds (see 
Paper I). We found that only the component at Visr =37 
km s _1 is sharply peaked in intensity at the IRAS posi- 
tion (19 h 31 m 12 s .8, 19°50' 22", B1950). Mapping of line 
emissions of high density tracers, HCN, HCO + , and CS, 
resulted in that the component at V\ ST =35-37 km s _1 is 
sharply peaked also at the IRAS position of this source. 
In this paper, we only show a result of the 5-point l'-grid 
cross mapping in the HCO + J = 1—0 transition in figure 
2. We can recognize from this figure that the H 12 CO + 
J= 1-0 line consist of the two components: a narrow com- 
ponent at V[ sr =37 km s _1 , and the weak broad component 
spreading from Vj sr =18 km s^ 1 to 50 km s~ x . Because 
these two components are also found in the other molec- 
ular lines as CO, HCN, and SO, and because they are 
strongly concentrated on the IRAS position of this object, 
we conclude that they are not from the fore/background 
clouds. Because strong lines as the CO J =1-0, 13 CO 
J =1-0, C 18 J = 1-0, CS J = 2-1, and HCN J = 1-0 
(Paper I), are somewhat contaminated by emission from 
the back/foreground molecular clouds, we presented the 
result of the mapping by these lines separately in Paper I. 
The results for HCO + described in this paper, i.e., the nar- 
row and broad components and their spatial extensions, 
were completely confirmed by mm-wave intcrfcrometric 
observations using the BIMA array (Nakashima, Deguchi 
2004) 

2.1. Time variation of maser lines 

We also observed maser lines to check their time varia- 
tion. The results are summarized in table 2 and the line 
profiles are shown in figure 3. It was striking that the SiO 
J = 1 v — 2 spectrum in 2001 February-March (mid- 
dle two panels in figure 3) exhibited strong emission at 
V\ SI ~ 50 km s _1 , but almost none in the J = 1-0 v = 1 
and J = 2-1 v = 1 spectra. The integrated intensity of 
the SiO J = 2-1 v = 1 line increased by a factor of more 
than 3 in about 40 days in 2001 March. If we compare 
these line profiles with those taken in May 2000 [figure 2 
in Nakashima, Deguchi (2000)], we can recognize drastic 
changes in the line profiles during 8 months. In May 2000, 
the radial velocities of most of maser lines were crowded 
around V\ SI ~ 20 km s~ x . However, we observed in 2001 
February that they spread between Vj sr ~ 20 km s _1 and 
55 km s _1 . This suggests that the center velocity of the 
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star is ~ 37 km s _1 , when we assume these maser lines are 
formed in the front and rear parts of the expanding shell. 

Top two panels of figure 3 show a drastic time vari- 
ation in H2O maser profiles during one- year time span. 
In 2001 February, emission appeared around the stellar 
velocity, ~ 37 km s _1 . While, after one year later, the 
profile looked doubly peaked and the line intensity be- 
came larger . Such a rapid change of the H2O maser line 
profile between the groups A (single peak at stellar veloc- 
ity) and group B (double peaks) in water masers (Engels 
et al. 1986) indicates that the direction of maser ampli- 
fication changes in a short time (Dcguchi 1977). If we 
assume that the double peaks of H2O masers are emit- 
ted in the front and rear part of the expanding envelope, 
we can derive V\ SI =33.6 km s _1 for the systemic veloc- 
ity and 17.3 km s _1 for the expansion velocity of the H2O 
maser shell. The intensity of the lower velocity component 
tends to surpass that of the higher velocity component in 
circumstellar H2O masers, indicating a maser amplifica- 
tion of continuum radiation of the central star (Takaba et 
al. 2001). 

Emissions of OH 1612, 1665, 1667 MHz lines have been 
detected toward this object (Lewis 2000, private commu- 
nication). Though the profiles of these OH lines are not 
exactly doubly peaked, emissions spread between Vi sr =10 
and 55 km s _1 . This observation also suggests that the 
stellar velocity is around Vj sr — 33 km s" 1 , which agrees 
with the stellar velocity derived from SiO and H2O profiles 
within uncertainty. 

2.2. Thermal lines of oxygen-bearing molecules 

From the detections of SiO, H2O and OH maser lines, 
we can safely suggest that the circumstellar envelope of 
this object is O-rich. However, the detections of H 13 CN 
J = 1-0 line toward this object and the velocity shift to 
the maser lines (Nakashima, Deguchi 2000) puzzled us, 
and these facts yielded the argument if the thermal lines 
may be due to contamination by background. 

To clarify the situation, we observed several oxygen- 
bearing molecules as SiO (J — 1-0, v = 0), SO, SO2, and 
HCO + . The results are summarized in table 3, and the 
line profiles are shown in figure 4. In addition, several 
lines of the various O-bearing molecules were searched 
for; frequencies of some transitions of these molecules were 
fallen in our spectrometer band coverages. We carefully 
checked whether or not any emission was seen in the spec- 
tra. These negative results are summarized in table 4. 

The line profiles of these O-bcaring molecules (except 
H 13 CO + ) are composed of two components: narrow sharp 
feature peaked at Visr ~ 36 km s _1 , and the broad feature 
with a full width at zero maximum of about 30 km s _1 
(see figures 2 and 4). It seems that the SiO J = 2-1 v = 
line profile does not have any narrow component. For 
the case of the S0 2 10(1,9)-10(0,10) line, the detection 
is marginal and it is difficult to tell whether it has the 
narrow or broad component. 

It is well known that the profile of the optically thin line 
in the expanding shell is flat-topped and that of the thick 
line is parabolic when the telescope beam cannot resolve 



Fig. 2. Line profiles of HCO+ and H 13 CO+ J = 1-0 tran- 
sitions. The observed position for each spectrum is indi- 
cated in the parentheses [ (AR.A., ADec.) relative to the 
IRAS position] on the right. The HCO + profile at the cen- 
ter clearly shows broad (V5 Br = 20-50 km s — 1 ) and narrow 
(around Vi sr = 36 km s _1 ) components. 
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Fig. 3. Time variation of the H2O and SiO maser lines. Dates of observations are shown in yymmdd format under the object name 
in each panel. 



the envelope (Morris & Alcock 1977). The profiles of 
the broad components in figure 4 do not seem to be flat- 
topped, but rather they look like parabolic. Therefore, we 
fitted the profiles of the broad component by a parabolic 
line-shape function (excluding the narrow part of emis- 
sion) and obtained the line parameters; the results are 
given in the 7 9 columns in table 3. The line parameters 
for the narrow components, which were computed from 
the residue by the parabolic fit, are also given in the last 
two columns in table 3. 

The average center velocity of the broad components, 
34.3 km s~ x , is very slightly shifted from the average peak 
velocity of the narrow component, 35.6 km s _1 . In addi- 
tion the profile of the SiO thermal line (J = 2-1 v =0) in 
figure 4 looks asymmetric with respect to the center veloc- 
ity. The difference in the radial velocity between the nar- 
row and board components seems significant. However, 
because the broad component is quite weak and over- 
lapped with the narrow component, it is difficult to sep- 
arate the broad component feature from the narrow com- 
ponent clearly in the present observations. Moreover, the 
profile of the CO J = 1-0 line (top right panel of figure 1 



of Paper I) clearly indicates a presence of absorptions on 
both high and low velocity sides of the narrow component, 
but with a slightly stronger absorption in the lower veloc- 
ity side. Therefore, the difference of the center-velocity 
between the narrow and broad components could be due 
to radiation transfer effect (absorption at lower velocity 
side) by a weakly expanding part of the cloud with low 
excitation temperature (for example, Morris et al. 1985). 

2.3. Thermal lines of C- and N -hearing molecules 

We have also detected a number of lines of C- and N- 
bcaring molecules. The results are summarized in table 5, 
and line profiles are shown in figures 5 and 6. These lines 
often have hyperfine splitting except the CS lines. The 
spectra of these lines do not seem to have the broad com- 
ponent but they exhibit the narrow component, although 
the hyperfine splitting (shown by arrows in figure 6) tends 
to conceal the presence of the broad component. 

Line widths of narrow components at the half intensity 
maxima are approximately 1.5-2.0 km s . The average 
center velocity of these lines in table 5 is 35.9 km s _1 , 
which is consistent with the average velocity of the narrow 
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Table 2. Time variations of maser line intensities. 
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Mol. Frequency Transition Vi sr (Peak) T*(Peak) Integ.Inten. rms Obs.date 
(GHz) (km s" 1 ) (K) K (km s' 1 ) (K) (yymmdd) 



H 2 


22. 


235 


6(16)-5(23) 


35 


.6 


2. 


,588 


12 


.105 





.031 


010210 


H 2 


22 


.235 


6(16)-5(23) 


1G 


.3 


6 


,628 


12 


.526 





.323 


021216 


H 2 


22 


.235 


6(16)-5(23) 


50 


.8 





,372 





.607 





.323 


021216 


SiO 


43 


,122 


1-0 v =1 


41 


.1 





,108 





.122 





.021 


010209 




42 


.821 


1-0 v =2 


49 


.9 





,184 





.817 





.021 


010209 




86 


243 


2-1 v =1 


36 


.9 





184 





.513 





.020 


010209 


SiO 


43 


,122 


1-0 v =1 

















.033 


010319 




42 


,821 


1-0 v =2 


53 


.8 





328 


1 


.347 





.042 


010319 




86, 


243 


2-1 v =1 


23 


.2 





199 


1 


.731 





.020 


010319 


SiO 


43 


,122 


1-0 v =1 


51 


.8 





213 





.376 





.055 


020519 




42 


,821 


1-0 v =2 


51 


.0 





,414 


1 


.036 





.067 


020519 


SiO 


43 


,122 


1-0 v =1 


50 


.8 





,186 





.346 





.030 


021215 




42, 


,821 


1-0 v =2 


51 


.7 





,420 





.970 





.043 


021215 



Table 3. Line parameters for line with broad and narrow features 



peak feature broad compo. narrow compo. 



Mol. 


Freq. 


Transition 


Hsr(p) 




rms 


Visr(par) T*(par) 


FWZM 


T*{nar 


) I.I.(nar) 




(GHz) 




(k„ B -l) 


(K) 


(K) 


(km a" 1 ) 


(K) 


(km a" 1 ) 


(K) 


(K km s -1 ) 


SiO 


86.847 


2-0 v=0 


33.3 


0.128 


0.027 


35.1 


0.086 


34.9 






H 13 CO+ 


86.754 


1-0 


36.2 


0.123 


0.021 








0.123 


0.556 


HCO+ 


89.189 


1-0 


36.9 


0.326 


0.025 


32.3 


0.060 


41.8 


0.269 


0.849 


SO 


99.300 


3(2)-2(l) 


35.8 


0.624 


0.019 


33.5 


0.112 


44.5 


0.513 


2.459 




109.252 


2(3)-l(2) 


35.6 


0.104 


0.043 


36.2 


0.028 


35.2 


0.076 


0.093 




100.030 


4(5)-4(4) 






0.023 












so 2 


104.239 10(1,9)-10(0,10) 35.8 


0.040 


0.012 













Table 4. Negative results for S-bearing molecules 



Mol. 


Freq. 
(GHz) 


Transition 


rms 
(K) 


so 2 


86.153 
86.639 
90.548 
104.029 


39(9,31)-40(8,32) 
8(3,5)-9(2,8) 

25(3,23)-24(4,20) 
3(l,3)-2(0,2) 


0.016 
0.025 
0.049 
0.020 


34 S0 2 


88.720 
104.391 


7(3,5)-8(2,9) 
10(1,9)-10(0,10) 


0.042 
0.025 


34 SO 


97.715 


3(2)-2(l) 


0.063 


ocs 


109.463 


9-8 


0.023 


SiS 


90.771 


5-4 


0.058 


C 33 S 


48.585 


1-0 


0.033 


C 34 S 


48.206 


1-0 


0.027 



component of O-bearing molecules. The width of the nar- 
row component seems to vary from molecule to molecule. 
The CS line has a full width of about 8 km s _1 , while, 
HCO + less than 2 km s _1 . This fact can be interpreted 
by a characteristic of the lines formed in a turbulent me- 
dia where the width depends on the line optical depth [for 
example, see Park & Hong (1995)]. 

The molecular species exhibiting the narrow compo- 
nent are somewhat similar to those observed in the C-rich 
circumstellar envelopes (e.g., Bujarrabal et al. 1994). A 
number of carbon-chain molecules have been detected in 
the C-rich protoplanetary nebulae (Fukasaku et al. 1994). 
However, SiS, which is common in the C-rich envelopes, 
was not detected in IRAS 19312+1950. 

The high HNC/HCN abundance ratio is often consid- 
ered as a signature of photochemical reactions (Morris et 
al. 1987; Glassgold 1996). The line intensity of the HNC 
emission, T* ~ 0.6 K, is comparable with the HCN inten- 
sity, T* ~ 0.6 K in this object, suggesting a photochemical 
origin of this molecule. Even taking into account the line 
saturation in the HCN lines [t(F = 2—1) ~2-4 can be 
deduced from T*(F = - l)/T*[F = 2 - 1) = 0.56 and 
T a( F = 1 - !)/ T a ( F = 2 - 1) = 0.76 in this star, where 
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Table 5. Line parameters for species with narrow features only 



Mol. 


Freq. 
(GHz) 


Transition 


Vl sr (peak) T*(pcak) 
(km s" 1 ) (K) 


Integ.Int. 
(K km s" 1 ) 


rms 
(K) 


NH 3 


23.694 


(1,1) 


35.1 


0.091 


0.731 


0.012 




23.723 


(2,2) 


37.4 


0.041 


0.203 


0.010 




23.870 


(3,3) 








0.010 




24.139 


(4,4) 


— 


— 


— 


0.022 


CS 


48.991 


1-0 


36.2 


0.621 


1.977 


0.042 




97.981 


2-1 


36.3 


1.100 


3.852 


0.059 




Oft A 1 Q 


O 1 

z— 1 


35.8 


0.191 


n A a 1 
U.441 


m 1 q 
U.Ullo 


H 2 CS 


104.617 


3(1,2)-(2(1,1) 


35.1 


0.046 


0.313 


0.010 


HNC 


90.600 


1 


35.4 


0.581 


1.930 


0.039 


N 2 H+ 


93.174 


1-0 


34.1 


0.122 


0.657 


0.020 




100.076 


11-10 


36.3 


0.163 




U.Ulo 


HC3N 


109.174 


12-11 


36.2 


0.104 


0.403 


0.020 


CN 


113.491 


1-0 (all) 


36.7 


0.521 


3.439 


0.033 


CH3OH 


96.739 


2(-l)-l(l) E 


35.4 


0.142 


0.384 


0.018 


CH3OH 


96.741 


2(0)-l(0) A+ 


35.9 


0.159 


0.448 


0.018 


c 17 o 


112.359 1-0 (+hyperfine) 


35.5 


0.197 


1.101 


0.027 


13 C 18 


104.711 


1-0 


36.4 


0.045 


0.200 


0.010 



these ratios are 0.2 and 0.6, respectively, in a thermal 
equilibrium], we can reach the same conclusion that the 
HNC/HCN abundance ratio is relatively high (> 0.3) in 
the narrow component. 

The molecular species exhibiting the narrow component 
are also similar to those found in the dark clouds; NH 3 , 
SO, H 2 CS , and CH 3 OH, etc., were observed in LDN 134N 
[sec Dickens et al. (2000)]. Therefore, in order to check 
the similarity in molecular abundances to dark clouds, we 
also searched for H2CS and CH3OH. We found emission 
from these molecules and the line profiles are shown in the 
bottom panels of figures 5 and 6. 

3. Discussion 

Paper I estimated the distance to this object as 2.5 kpc; 
the bolometric correction for the IRAS 12 micron flux 
density, {BC)\2 =3.8 (Deguchi et al. (2001)), was applied 
from the R (DSS), J, H, and K (2MASS) magnitudes and 
IRAS 12, 25, 60, and 100 fim flux densities, and typical 
luminosity of the AGB object, 8000 Lq, was used. For 
simplicity, we adopt this value for a distance of this object 
in this paper. 

Paper I gave two reasons for the central star to be an 
evolved object: (1) the presence of deep CO first overtone 
absorption bands at 2.3 (im as usually seen in late Mill 
stars, and (2) no nearby starforming activity such as of- 
ten seen in YSOs. Furthermore, we can add the OH maser 
line property of this object: the dominant OH 1612 MHz 
line accompanying weak OH 1665/1667 HMz lines (Lewis 
2001; private communication), which is often used to sep- 
arate evolved objects from YSOs (Caswell 2000). All of 



these facts strongly support the evolved star interpreta- 
tion of the central star. 

There is evidence for O-rich material surrounding a 
class of carbon stars (Willems, de Jong 1986; Nakada 
et al. 1987; Le Bertre et al. 1990; Jiang 2000). These ma- 
terials arc considered to be the gas and dust which were 
ejected in the O-rich AGB era of the central star, but 
were trapped by a hypothetical binary companion (Lloyd- 
Evans 1990). When the central star changes the atmo- 
spheric composition from O-rich to C-rich, the presence of 
the O-rich material would observationally become evident 
due to their chemistry difference. In fact, a sharp narrow 
CO line profile (~ 1 km s _1 ) having slightly broader wing 
emission (~ 5 km s _1 ), which is similar to the narrow and 
broad components of IRAS 19312+1950, has been found 
in the silicate carbon star, BM Gem (Kahane et al. 1998; 
Jura, Kahane 1999). In order to explain O-rich mate- 
rials surrounding silicate carbon stars, a model of orbit- 
ing molecular reservoirs has been proposed (Lloyd-Evans 
1990; Jura, Kahane 1999). In contrast, Kerschbaum et 
al. (1999) found O-rich examples (semircgular and irreg- 
ular variables) of the narrow (< 5 km s" 1 ) feature su- 
perposed on the broad (~ 20 km s _1 ) component in CO 
lines. These may be interpreted as a result of cither mul- 
tiple winds or a Keplerian disk (Bergman et al. 2000). 

It is possible that we are looking these trapped mate- 
rials as a narrow-component cloud [for example, Kahane 
et al. (1998)]. Because the narrow and broad components 
found in the IRAS 19312+1950 line profiles have a simi- 
larity with the CO line features found in these stars, the 
model of orbiting molecular reservoirs must be pursued 
first in the present case too. However, as discussed in 



No. ] Circumstellar Chemistry of the Bipolar Nebula IRAS 19312+1950 7 



Table 6. Negative results for miscellaneous species 



Molecule Frequency Transition 
(GHz) 


rms 
(K) 


Freq. Ref. 


MgNC 


23.859 


2-1- 


0.006 


Kawaguchi et al. 1993 




23.875 


2-1+ 


0.006 


Kawaguchi et al. 1993 




47.726 


4-3- 


0.023 


Kawaguchi et al. 1993 




47.741 


4-3+ 


0.023 


Kawaguchi et al. 1993 




95.454 


N = 8—7- 


0.016 


Guelin et al. 1995 




95.469 


N = 8—7+ 


0.016 


Guelin et al. 1995 


26 A1C1 


89.351 


6-5 


0.016 


calc. from Guelin et al. (1995) 




104.237 


7-6 


0.013 


calc. from Guelm et al. (1995) 


CH 3 OH 


48.372 


l(0)-0(0) A+ 


0.028 


Lovas 1992 




48.376 


l(0)-0(0) A+ 


0.028 


Lovas 1992 


C 4 H 


104.667 


11-10 J=23/2-21/2 0.010 


Lovas 1992 




104.705 


11-10 J=21/2-19/2 0.010 


Lovas 1992 


H 2 Q 


96.261 


4(4,0)-5(3,3) u 2 


= 1 0.018 


Lovas 1992 



f ig4 . eps 



Fig. 4. Line profiles of SiO J = 2-1 v =0, SO 3(2)-2(l) and 
2(3)-l(2), and SO2 10(1,9)-10(0,10) transitions. The SiO 
profile exhibits only a broad component, and SO profiles ex- 
hibit both broad and narrow components. Because of the 
noise, it is hard to judge whether it is a narrow or broad 
component for SO2. 



the following sections, the derived total mass of the gas 
(~ 25 Mq within 15" radius), which is necessary to ex- 
plain the narrow component of IRAS 19312+1950, seems 
too large as an orbiting molecular reservoir. 

3.1. The broad component: molecular abundances 

In order to obtain molecular abundances, numerical 
calculations based on the large-velocity-gradicnt (LVG) 
model were performed. Here the broad-component emis- 
sion was assumed to come from the expanding envelope 
of an AGB/post-AGB object. The numerical code devel- 
oped by Dcguchi et al. (1990) was used to compute line 
intensities of linear molecules in the expanding shell. The 
code was modified to involve SO molecule which has a 
slightly complex energy-level structure; transition prob- 
abilities (Einstein A coefficients etc.) of SO were taken 
from Tiemann (1974) [see the level diagram in Omont et 
al. (1993)]. In this numerical code, the kinetic tempera- 
ture of the envelope is calculated by taking into account 
the adiabatic and line (CO, CI, and CII) cooling and the 
gas-dust drift and CI ionization heating (see the detail 
in Deguchi et al. 1990). Dissociations of CO, HCN, SiO, 
and SO due to interstellar UV radiation are taken into ac- 
count; self absorptions of CO UV lines are also taken into 
account. However, photodissociation rates of SiO and SO 
are not well known. Therefore, we assumed that the pho- 
todissociation rates are the same as the CO dissociation 
rate (without self-absorption). 

We assumed that the envelope is expanding with a 
constant velocity of 25 km s -1 . The calculation starts 
from the inner radius of 1.5 x 10 16 cm and stops at the 
outer radius of 1.87 x 10 17 cm (corresponding to a 5" ra- 
dius from the central star). The mass loss rate (to be 
constant), kinetic temperature at the inner radius, and 
molecular abundance were free parameters and the best- 
fit model to the observational data was sought. The re- 
sults are given in table 7, where molecular species, abun- 
dance of the observed molecule with the broad compo- 
nent, frequencies, transitions, and the observed main- 
beam brightness temperature, and the calculated main- 
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Fig. 5. Line profiles of HC 3 N, CS and H 2 CS. The hypcrfine 
splitting of HC3N is quite small and not recognizable here. 
Both of HC3N and CS lines seem to exhibit a narrow compo- 
nent. 



f ig6 . eps 



Fig. 6. Line profiles of NH 3 , HNC, N 2 H+, and CH3OH. The 
positions of the hyperfine lines were indicated by arrows. The 
splitting of HNC are less than 0.7 km s _1 , so that they are 
not shown here. All of these lines exhibits only a narrow 
components. 
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Table 7. Model fit for the broad component for mass loss 
rate of 2.6 X lCr 4 A/ yr" 1 



Mol. 


Abundance Frequency Transition 
(per H 2 ) (GHz) 


MB 

(K) 


T^modcl 
MB 

(K) 


SiO 


1.1 x 10-' 


86.847 2-1 v = 


0.18 


0.18 


HCN 


3.9 x 10~ 7 


88.632 1-0 


0.41 


0.41 


HCO+ 


5.7 x 1CT 8 


89.188 1-0 


0.12 


0.12 


"CO- 


4.0 x 10~ 6 


110.201 1-0 


<0.1 


0.09 


CO 


1.3 x 10~ 4 


115.271 1-0 


1.79 


1.79 


so 


5.5 x 10"' 


86.094 2(2)-l(l) 


0.05 


0.04 






99.299 3(2)-2(l) 


0.22 


0.20 






109.252 2(3)-l(2) 


0.05 


0.10 



f ig7 . eps 



Fig. 7. Line profiles of NH 3 , C 34 S, C 17 0, and 13 C ls O lines. 

beam brightness temperature were shown for the model 
with the mass loss rate of 2.6 x 10 _4 Mq yr _1 . In this 
model, the main-beam brightness temperature was deter- 
mined mostly from line brightness temperature at around 
10 17 cm with Nh 2 = 1.1 x 10 4 cm" 3 (except the CO J—l- 
transition). Because photodissociation of molecules oc- 
curs slightly outside of this radius, a slight uncertainty 
of the photodissociation rate in various molecules except 
CO does not influence much to the results. The obtained 
SO abundance, 5.5 x 10~ 7 , is similar to those found in 
normal O-rich envelopes (Guilloteau et al. 1986;Omont et 
al. 1993). 

For the case of HCO + , the abundance of this molecule 
is assumed to be constant through the envelope in the 
present calculation. Because the molecular abundances 
of HCN and HCO + in the O-rich circumstellar envelopes 
were somewhat controversial (Morris et al. 1987; Deguchi 
et al. 1990), we also computed the expected HCO + abun- 
dance in the same code used in Deguchi et al. (1990). The 
maximum HCO + abundance can reach to 1.5 x 10 per 
H2 at the radius 2.3 x 10 17 cm (if mass loss occurs con- 
tinuously beyond 1.87 x 10 17 cm). Because there are still 
some arguments in the H^" dissociative recombination rate 
(Amano 1988), this value is probably available maximum. 
The obtained HCO + abundance is less than the value ob- 



tained from the present photo-chemical model (see also 
Glassgold 1996). Furthermore, dumpiness in the outflow- 
ing envelope strongly influences the photo-chemical prod- 
ucts (Redman et a. 2003). The spherical model with outer 
radius of 1.87 x 10 17 cm may be inadequate for explaining 
HCO + . In fact, a bipolar structure has been found for the 
HCO + emission distribution (Nakashima, Deguchi 2004). 

3.2. The narrow- component cloud: molecular abundance 

In order to obtain molecular abundances in the nar- 
row component cloud, we performed numerical calcula- 
tions under the LVG model. Though the LVG model was 
a choice for a matter of convenience here, more sophisti- 
cated, but formidable radiative-transfer model involving 
dumpiness and turbulence with velocity gradient (Park 
& Hong 1995) must be pursued for the narrow compo- 
nent cloud as a next step; there is evidence that the spec- 
tral line of the narrow component is broadened by turbu- 
lence as mentioned in section 2.3. However, the present 
data, which was obtained in a single-dish observation, do 
not provide much information on spatial distribution of 
molecular emission. Then, it seems meaningless to ar- 
gue on differences of outcomes between the LVG and tur- 
bulence models. Therefore, we used the LVG model for 
the narrow component cloud with inner and outer radii, 
1.87 x 10 17 cm and 5.6 x 10 17 cm (corresponding to 5 
and 15" radii), respectively The uniform expansion rate, 
dV/dr = 3.2 x 10 -13 s _1 , was adopted giving an expansion 
velocity of 2 km s _1 at the outer boundary. Radiative 
interaction with the central expanding core (broad com- 
ponent) was also taken into account as if the core emits 
continuum emission, but it turned out to be negligible. 
The number density was assumed to decrease with ra- 
dius as r -3 , and the absolute value of the density and the 
molecular abundances were adjusted to fit calculated to 
observed line intensities. For simplicity, the calculations 
were restricted for linear- type molecules (except SO). We 
solved level populations up to J = 6 for CO, HCN, HCO + , 
and HNC, to J = 8 for CS , and to J = 15 for HC 3 N at 
different radii and calculated the brightness temperatures. 
Multiplying the telescope beam pattern, we obtained the 
molecular abundances to fit the calculated to the main- 
beam brightness temperatures of observed lines. Kinetic 
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Table 8. Molecular abundance for narrow component 



Molecule Abundance Freq. Transition 





(per H 2 ) 




(GHz) 


(K) 


(K) 




CO 


1.3 x 10" 


— i — 

-4 


115.271 


1-0 


6.83 


6.83 


HCN 


4.0 x 10" 


- 9 


88.631 


1-0 


0.77 


0.77 


HCO+ 


1.2 x 10" 


. f) 


89.189 


1-0 


0.53 


0.53 


HNC 


3.4 x 10" 


-9 


90.663 


1-0 


1.12 


1.12 


HC 3 N 


7.0 x 10" 


-9 


100.076 


11-10 


0.31 


0.31 








109.174 


12-11 


0.19 


0.22 


cs 


1.0 x 10" 


-8 


48.991 


1-0 


0.85 


1.23 








97.981 


2-1 


1.77 


1.66 


so 


5.5 x 10" 


-9 


86.094 


2(2)-l(l) 


0.05 


0.27 








99.300 


3(2)-2(l) 


0.77 


0.77 








109.252 


2(3)-l(2) 


0.17 


0.42 



temperature in the model (uniform) was determined to 
be 13 K by adjusting the calculated 12 CO line intensity 
to the observed intensity. The results are shown in table 
8. It suggests that the abundances of HCN and SO are 
two orders of magnitude lower than those for the broad 
component (table 7), making the discrete difference of the 
molecular composition at the boundary between narrow 
and broad components. The total mass of the narrow 
component cloud in this model is 24.4 M© (within 15" 
radius), which is consistent with the previous estimation 
of the mass of this dust cloud based on a 18 CO column 
density (Nakashima ct al. 2004). 

3.3. isotopic abundance 

Another important clue to guess the origin of cool mate- 
rial surrounding the object is isotopic abundance. For this 
purpose, we obtained the isotope abundance ratios for the 
narrow-component cloud with the same LVG model used 
in the previous section. Table 9 summarizes the isotope 
line-intensity (peak-intensity) ratios and calculated abun- 
dance ratios (in the 2nd and 5th columns, respectively). 
For comparison, terrestrial and interstellar isotope abun- 
dance ratios are given in the last two columns. If lines are 
optically thin, the line intensity ratio approximately gives 
the isotope abundance ratio. However, because strong 
lines, as 13 CO J = 1-0, HCO+ J = 1-0 and CS J = 2-1, 
tend to be saturated in intensity (optically thick), they 
only give a lower limit of abundance ratio. 

The ratio, C 18 0/ 13 C 18 0, is very often regarded to give 
the 12 C/ 13 C abundance ratio in dark clouds in a good ap- 
proximation [for example, sec Langcr, Penzias (1993)] 
because these lines are optically thin. In fact, in the 
LVG model computations, the line intensities of the mi- 
nor CO isotopic species vary almost proportional to the 
abundance of CO [times (dV/dr) -1 too] [e.g., figure 17 of 
Castets et al. (1990)], as far as it is optically thin. We ob- 
tained the ratio 12 C/ 13 C~ 27 ±9 from the C 18 0/ 13 C 18 
line ratio in table 9. It is well known that 13 C is enriched in 
envelopes of some carbon stars [for example, see Woods et 
al. (2003)]. Furthermore, Jura et al. (1988) found low line 



intensity ratio of 12 CO to 13 CO J =1-0 transition in sev- 
eral optical carbon stars which exhibit 13 C-richness. The 
observed abundance ratios are approximately between 3 
and 80 for evolved objects [for example, Schoer, Olofsson 
(2000)]. No enrichment of 13 C compared with the case 
of interstellar clouds (in table 9) indicates that the mate- 
rial is unlikely to be an ejecta from a (supposedly) dead, 
C-rich AGB companion. 

However, we have to remind that the detection of 
13 C 18 J = 1-0 line was marginal (see the lowest panel 
of figure 7). Taking into account the uncertainty of the 
temperature calibration and pointing errors of the 45-m 
telescope at 110 GHz due to wind etc., we think that the 
observed value involves uncertainty of a factor of about 2. 
Furthermore, if we consider the model dependency of the 
CO line intensity, especially in turbulent clumpy clouds, 
and of emission size difference in C 18 and 13 C 18 lines, 
we think that the conclusion is somewhat provisional. 

Kahane et al. (1992) observed several carbon stars and 
planetary nebulae, and found that the 18 0/ 17 ratio is 
less than unity. However, the present case is apparently 
different from such C-rich objects; the C 18 0/C 17 ~ 5 in 
table 9. Concerning on 17 0, we can safely conclude that 
it is not appreciably enriched in the narrow-component 
cloud as in these C-rich objects. Rather, the 18 0/ 17 
ratio coincides with the interstellar value. In conclusion, 
we found no evidence of the isotope enrichment /deficiency 
compared with the interstellar values in the present ob- 
servation. 

3.4- Narrow component: what is it ? 

The detections of SO, HCN, HCO+, and NH 3 first 
suggested a resemblance of chemistry between IRAS 
19312+1950 and OH231. 8+4.2 (a proto-planetary neb- 
ula), leading a conclustion that IRAS 19312+1950 was 
a protoplanctary nebula. In fact, the first three 
molecules were observed to have a broad component as 
well. However, we recognized several differences between 
the molecular species found in IRAS 19312+1950 and 
OH231.8+4.2. Firstly, SO2 molecules are quite deficient 
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Table 9. Isotope line-intensity and abundance ratio (narrow component) 



Molecule 


Peak I. ratio Uncertainty Transition Calc. 


abu. 


ratio Terrestrial ratio 1 


Interstellar ratio 


i3 CO/ i3 C i8 


66 


± 20 


1-0 


498 


489 


489 2 


13 CO/C 18 


3.6 


± 0.3 


1-0 


19 


5.43 


7-25 3 


C 18 /13 C 18 


18 


± 6 


1-0 


27 


89.1 


20-70 4 


c 18 o/c 17 o 


4.1 


± 0.7 


1-0 


5.7 


5.5 


3.2 2 


HCO+/H 13 CO+ 


2.2 


± 0.5 


1-0 


10 


89.1 


~50 5 


cs/c 34 s 


4.8 


± 1.0 


2-1 


9.1 


22.5 


24 5 



1 from appendix VII of Townes, Shawlow (1975). 

2 Wilson, Matteucci (1992) 

3 Gredcl ct al. (1994) 

4 Langer, Penzias (1993), Savege et al. 2002; see references therein. 

5 Lucas (1998) 

6 Chin et al. (1996) 



in IRAS 19312+1950 (sec table 4), but were abundant 
in OH231. 8+4.2 (Morris et al. 1987), although SO lines 
are detected in both. The line intensity ratio of the SO2 
3(13)^2(02) (104.029 GHz) to SO 2(2)-l(l) line (86.094 
GHz) is 6.2 in OH231. 8+4.2 (Morris et al. 1987); the abun- 
dance ratio of S0 2 to SO is about 4-7 (Claude et al. 2000), 
and somewhat similar values have also been found in O- 
rich AGB stars (Sahai, Wannier 1992) (except 2 super- 
giants, NML Cyg and VY CMa, which seem to be S0 2 - 
poor relative to SO; Omont ct al. 1993). However, this 
intensity ratio is < 0.2 in IRAS 19312+1950 (in the condi- 
tion that the marginally detected SO2 line is regarded to 
be a narrow component). The deficiency of SO2 relative 
to SO in IRAS 19312+1950 could be attributed to a differ- 
ence in O/C abundance ratio between two stars, or to less 
developement of shock compared with OH231.8+4.2 (e.g., 
Claude et al. 2000). Secondly, the NH 3 (1,1) and (2,2) in- 
version lines were detected in IRAS 19312+1950. The 
line intensity ratio suggests the rotational temperature of 
the concerning levels (both para-NHa) of 19 K for the 
narrow component cloud. For the case of OH 231.8+4.2, 
the line intensities suggested the NH3 rotational temper- 
ature of (27 + 4) K (Morris et al. 1987), indicating that 
the narrow-component cloud is slightly cooler than the 
envelope of OH231. 8+4.2 

Rich molecular lines detected toward IRAS 19312+1950 
(as a narrow component) rather resemble to the lines usu- 
ally found in dark clouds. For example, in L134N, the 
SO lines have been found to be strong but not much SO2 
(Swade 1989; Dickens et al. 2000). Therefore, the de- 
ficiency of SO2 to SO can be explained naturally if we 
consider that the narrow-component cloud is a dark cloud 
as L134N. 

The detection of CH3OH in this paper also favors the 
dark-cloud interpretation of the narrow component cloud. 
The detected CH 3 OH lines at 96.74 GHz are thermal 
emission and have been found widely in dark (Friberg et 
al. 1988) and translucent clouds (Turner 1998). Note 
that a search for the 6.7 GHz CH3OH 5i-6o A + line, one 
of the strongest CH 3 OH maser lines, has been negative 
in IRAS 19312+1950 (MacLeod et al. 1998; Szymczak et 



al. 2000), although numerous maser lines of CH3OH 
have been detected normally in dense massive star forming 
clouds (in this sense, the cloud toward IRAS 19312+1950 
is not a massive star forming region). Though CH3OH 
could be produced in circumstellar envelopes of O-rich 
stars (Charnley et al. 1995), it has never been detected 
in evolved stars to the authors' knowledge [see Latter, 
Chanley (1996) ; Charnley, Latter (1997)]. 

Time-dependent chemical models of low-mass dark- 
cloud cores (Nejad, Wagcnblast 1999) gave a result that 
abundances of C-bearing complex molecules as CH3OH 
and CH3CN are peaked at a few times 10 5 years, while 
the peak abundance of SO and SO2 are attained at later 
times of ~ 10 6 years (though these time scales are depen- 
dent on the density of the cloud). Note that the chemistry 
of CH3OH and CH3CN in translucent clouds is still quite 
controversial [see an excellent summary in Turner (1998) 
on gas-phase and grain-surface chemistry], and the con- 
finement of their chemical evolution at relatively earlier 
time in the Nejad, Wagenblast (1999)'s model was criti- 
cally affected by cosmic-ray-induced photoionization and 
photodissociation processes. 

It is useful to compare the time scale of molecule forma- 
tion with the other physical time scales of the envelope. 
The time scale of expanding gas with 10 km s _1 crossing 
the envelope (15" ~ 5.8 x 10 17 cm) is as short as 1.8 x 10 4 
y. Therefore, even if we interpret the narrow-component 
cloud being the stellar ejecta, we need relatively long-lived 
(~ 10 6 years) trapped materials to create CH3OH. The 
period of circulation around a 5 M Q star at 10" radius is 
1.8 x 10 6 y. Therefore, it is possible to create the complex 
molecules as CH3OH within one orbital motion of ejected 
gas. Because the life time of the thermal-pulse AGB era 
is < 10 6 y for a progenitor mass of > 4 M Q (Blocker 1995), 
the ejecta must come from the companion, not from the 
AGB/post-AGB star which is seen now as a central star 
of the mass of about 5 Mq (Paper I). 

It is possible to consider that material of the narrow- 
component cloud is a remnant of the dark cloud which en- 
countered with this object several million years ago. Such 
encounter with a dark cloud is considered to occur rela- 
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tively frequently in the solar neighborhood (in every ~10 
Myr per a star; Griv et al. 1997; Fresneau et al. 2003). 
Taking into account the chemical-evolution time scale of 
dark clouds, we estimate that such encounter must have 
occurred in the last 10 6 yr ago for IRAS 19312+1950. In 
this case, the encountered cloud cannot be very far away 
from the object at present; for example, it must be lo- 
cated within 5' from IRAS 19312+1950, provided that 
the relative velocity between the star and the cloud is 10 
km s . This is also an attractive hypothesis in view of 
an explanation for silicate carbon stars and narrow CO 
line profiles of short-period variable stars. However, the 
large mass of the narrow-component material (as much as 
~ 20A/q) again causes a problem in this interpretation; It 
must be quite difficult to accrete the mass more than the 
star mass during the passage of such dark cloud and to 
keep it around the object during a few Myr. 

Note that brightness of the central star at near-infrared 
bands (K = 6.3 in IRAS 19312+1950 ; Paper I) is so dif- 
ferent from class-0 YSOs as L134N or L1551 IRS5. These 
objects usually do not have any bright near-infrared coun- 
terparts. For example, L1551 IRS5 does show very faint 
patchy structure in the K-band images [Hodapp 1994; or 
see Momose et al. (1998) for the overlay image with the 
C 18 J = 1-0 map] . As far as we see 2MASS images of the 
similar type of objects as L134N (Dickens et al. 2000) and 
IRAS 16293-2422 (Blake et al. 1994), we find no near- 
infrared counterpart. Though HH objects as B5 IRS 1 
(Fuller et al. 1991; K=11.21) tend to have brighter near- 
infrared counter parts (Reipurth, Bally 2001), these ob- 
jects have well developed CO outflow too. On the other 
hand, the size of high-velocity flow of IRAS 19312+1950 
(if the broad component is interpreted as a bipolar flow ob- 
served in YSOs) is small compared with the size of narrow 
component cloud (~ 20"), suggesting again a difficulty of 
such interpretation that the central star is a young stellar 
object. 

From above arguments, we are forced to conclude that 
the narrow-component cloud is a small dark cloudlet. 
However, because the presence of a ring structure and 
spurs seen in figure 1 indicates a physical relation between 
the central star and the surrounding cloud, the narrow 
component must be a part of the object, which is phys- 
ically associated with the central star, even though they 
may not be the gas ejecta from the central star or hypo- 
thetical companion. In order to check the possibility of 
YSO further, it might be nice to observe CO and CO2 ice 
absorption bands at 4.2-4.7 /im (as a revelation of molecu- 
lar clouds: Gibb et al. 2000) and crystalline H 2 at 44/62 
/im (as a revelation of post-AGB outflows; Maldoni et 
al. 2003) for the future work . 

It is probably the best to perform mm-wave interfero- 
metric observations with high spatial resolution in order to 
clarify the physical relation between the narrow and broad 
components. Preliminary results of intcrfcromctric map- 
ping of the CO J = 1-0/2-1 emissions with Nobeyama and 
Bcrkclcy-Illinois-Maryland Arrays gave that the broad 
component was partially resolved in ~ 3" beam, but the 
narrow component is extended to about more than 10" in 



size. These results are consistent with the 45-m mapping 
observations (Paper I), and they fit well to the present 
interpretation that the broad component is the expand- 
ing envelope of the AGB/post-AGB object which is sur- 
rounded with a cool cloud at the distance of 2.5 kpc. 

4. Conclusion 

We have detected a variety of molecular lines of 0-, 
C- and N-bearing molecules toward IRAS 19312+1950, a 
bipolar nebula. The line profiles consist of narrow and 
broad components of the widths of ~4 and ^30 km s _1 , 
respectively. From the parabolic line shape, the broad 
component is inferred to originate from the circumstcllar 
envelope of this object. The LVG calculations suggest 
the mass loss rate of about 2.6 x 10 _4 Mq yr -1 for this 
object. However, the origin of the narrow component is 
still a puzzle; from the low temperature and variety of 
molecules detected, it is suggested to be a dark cloudlet 
with a mass of about 20 M@. If it were a trapped material 
by a hypothetical companion, the H2/CO mass ratio must 
be a factor of 10 less than the value normally assumed. If 
this is really a dark cloud with a mass more than 20 Mq , 
how the AGB/post-AGB star can accommodate a cloud 
more massive than the star mass ? Furthermore, if this 
is a low-mass young stellar object with SiO masers more 
than a few kpc away, it is hard to explain the apparent 
size of the bright nebulosity and luminosity of the central 
star. 

It is hard to provide a reasonable interpretation for this 
object except considering somewhat exotic, previously un- 
known, star model as a hypergiant (Jura et al. 2001), or 
possibly an AGB/post-AGB object with the huge orbiting 
gas reservoir in which the H2/CO ratio is small compared 
with a normal ratio. Another extreme is a low-mass young 
stellar object with SiO masers with yet undeveloped CO 
outflow but shining the surrounding dust of the size of 
30". Neither of them could be acceptable, except with 
excuse if they were extreme objects. 

The authors thank Dr. M. Tamura, and K. Murakawa 
for providing them with the SIRIUS image of IRAS 
19312+1950. They also thank Dr. T. Umemoto for useful 
suggestions, and students (H. Fukushi, and Y. Tamura) 
for help of observations. 
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